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We study the coarsening of two-dimensional oblique stripe patterns by numerically solving po- 
tential and nonpotential anisotropic Swift-Hohenberg equations. Close to onset, all models exhibit 
isotropic coarsening with a single characteristic length scale growing in time as t^^^. Further from 
onset, the characteristic lengths along the preferred directions x and y grow with different exponents, 
close to 1/3 and 1/2, respectively. In this regime, one-dimensional dynamical scaling relations hold. 
We draw an analogy between this problem and Model A in a stationary, modulated external field. 
For deep quenches, nonpotential effects produce a complicated dislocation dynamics that can lead to 
either arrested or faster-than-power-law growth, depending on the model considered. In the arrested 
case, small isolated domains shrink down to a finite size and fail to disappear. A comparison with 
available experimental results of electroconvection in nematics is presented. 

PACS numbers: 47.54.-r, 47.55.P-, 61.72.Cc 



I. INTRODUCTION 

Domain coarsening of uniform phases in systems lo- 
cally in thermal equilibrium has received a lot of at- 
tention over the past decades [l|, After quenching 
a system below a transition point, ordered regions of 
the lower symmetry phase form and their characteris- 
tic length L grows slowly, usually as a power law at large 
times, L ^ t". Numbers of theoretical schemes based on 
the dynamical scaling hypothesis and on the study of the 
relaxation of topological defects have been proposed to 
infer coarsening exponents and address their universality 

Biiiia. 

Less understood is the phase ordering kinetics in sys- 
tems driven out-of-equilibrium, in particular those form- 
ing dissipative structures like regular patterns. Ba- 
sic non-equilibrium structures are the modulated phases 
(stripes) of well defined periodicity that appear above an 
instability threshold in thermal convection of fluids, in 
driven chemical reactors or in electroconvection of ne- 
matic liquid crystals [1, H, [13, [ll| . 

A body of numerical results [H, [H, [H, [H, M, M, 
H E, |2fl IM in m, [13, HI suggests that the order- 
ing kinetics of stripe patterns following a quench is not 
universal but model-dependent, and even, parameter- 
dependent. Modulated phases in isotropic systems ac- 
tually exhibit a rich variety of topological defects (dis- 
locations, disclinations, grain boundaries) not observed 
in uniform phases. There might exist many mechanisms 
controlling the coarsening of stripe patterns: for instance, 
the annihilation of disclination quadrupoles, as observed 
experimentally [2^,[23,[2l| in block-copolymer melts with 
cylindrical mesophases (leading to n — 1/4), or grain 
boundary motion in the Swift-Hohenberg model close to 
onset (n = 1/3) [13, [H US Hi ■ 
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For a given system, the growing length scale may even 
be not unique and the dynamical scaling hypothesis not 
fulfilled. Different defect-based definitions of the size L of 
the ordered striped domains can lead to different coars- 
ening exponents [3, [H, Hi], . This effect seems to be 
more pronounced far from the onset of stripe formation. 
In this regime, the growth of order often becomes very 
slow or arrested (15^] , in particular due to pattern- induced 
pinning effects acting on defects [13, HI, HI] • 

In addition, unlike equilibrium phases, the dynam- 
ics of nonequilibrium patterns is a priori not driven by 
the minimization of a free-energy functional but bet- 
ter modelled by nonpotential (nonrelaxational) equations 
d, [13, [O, [HI- Recent simulations of the Boussinesq 
equations for Rayleigh-Benard convection [l^l suggest 
that multiple length scales must be defined to char- 
acterize the coarsening of nonpotential stripe patterns. 
Adding nonpotential terms to potential models [13, [s^l 
(or, similarly, taking into account hydrodynamics in the 
order parameter description of block-copolymer melts 
[25| ) can noticeably increase some coarsening rates and 
lead to qualitative changes in the structure of the de- 
fected patterns. 

In this context, systems forming anisotropic patterns 
have recently attracted a particular attention [3l|, [S^, [23, 
[33,[3H|- An example are the obhque phases, composed of 
stripes that make a fixed angle, 9 or —6, with respect 
to some X-axis. Domain coarsening in these systems 
have been studied in details experimentally in electro- 
convection of nematics jsil . [s^ . [33j . Anisotropic stripes 
may be regarded as topologically simpler than isotropic 
ones because of their finite number of orientations. How- 
ever, some challenging difficulties mentioned above for 
isotropic systems still remain, such as the presence of 
various sorts of defects and nonpotential effects. 

The experimental studies of refs. [3l|, [s^l on oblique 
stripes report that chevron grain boundaries and dislo- 
cations arrays (the domain interfaces along the x and 
y directions, respectively) have very different mobilities. 
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After a quench, their respective densities have different 
power-law decays with time, both with unusually small 
exponents. These results remain largely unexplained the- 
oretically. A first numerical study '34!| based on a po- 
tential anisotropic Swift-Hohenberg equation for oblique 
stripes jsH] reproduced the decay of the dislocation den- 
sity but predicted a much faster ordering in the y direc- 
tion than observed experimentally. 

We revisit this problem by introducing in Section II 
nonpotential versions (never considered before, to our 
knowledge) of the anisotropic model used in [3J]. The 
wavenumbers selected in these models are comparable to 
the experimental ones. In Section III, we define an ori- 
entational order parameter, its correlation lengths along 
the preferred x and y directions, and the defect densities 
(chevron and dislocations). In Section IV, we find two 
robust coarsening regimes at moderate quench depths, in- 
dependent of the potential or nonpotential nature of the 
models: an isotropic regime very close to onset and an 
anisotropic one (in partial agreement with experiments) 
at intermediate quenches. One-dimensional dynamical 
scaling relations hold in the latter case. In section V, we 
show that, up to moderate quenches, the problem can 
be approximately reduced to a Model A in a stationary, 
spatially modulated external field. Far from onset, non- 
potential effects can lead to arrested or very fast domain 
growth, depending on the model. In the arrested case, 
small domains only shrink down to a finite size. Conclu- 
sions are presented in Section VI. 



II. BASIC MODELS AND THEIR FIRST 
ORDER AMPLITUDE EQUATIONS 

An electroconvection set-up consists of a doped ne- 
matic liquid crystal confined between two plates with pla- 
nar alignment. When an ac transverse electric field is ap- 
plied, normal or oblique stripes (among other patterns) 
can form depending on the frequency [13, [H, [H, [ioj . 
Some time ago, Pesch and Kramer (PK) introduced a 
phenomenological Swift-Hohcnberg-like model that ex- 
hibits a transition from normal to oblique rolls [36j : 

d,^ = rV--C'(V2 + fc2)2V;-^a>+^9X'^ + iVLM, 

(1) 

where '0(r, t) is the local order parameter, NL[ip] = 
c and rj are dimensionless anisotropy parameters, fco a 
characteristic wavenumber, C, (set to I/Zcq in the follow- 
ing) the coherence length and r the main control param- 
eter. In electroconvection of nematics with planar align- 
ment, X is the coordinate in the direction parallel to the 
undistorted director. 

Linear stability analysis of the uniform state ?A(r', t) = 
of Eq.(IT]) against small periodic perturbations 
dijj{r,t) = (5?/'oexp(ifc • r + at) gives the dispersion re- 



lation 

where k = px + qy. Maximizing ^ with respect to p 
and q and looking for the values of r for which ct = 0, 
one finds that some modes of finite wavenumber become 
marginally unstable when the control parameter r ex- 
ceeds the critical values ri"^ or ri°' : 

normal stripes: r'"'' = 0, {p^ = k'^,q^ ~ 0}, (3) 

oblique stripes: r^°^ = j+^i^ < 0' {Pc = ^+i|^^o. 

'I'c - TT^kl}^ (4) 

where we have assumed that c > 0. The marginally 
unstable wavevectors in the oblique case make an 
angle 6 or —9, with respect to the a;-axis, where 
9 = arctan[-y/?7/(c -t- rf)]. These two-degenerate oblique 
modes, noted as fc+ (zig) and fcj (zag) in the following, 
are observed for 77 > only. The model equation ^ ac- 
tually reproduces the transition from normal to oblique 
stripes when 77 is tuned from negative to positive values. 

In the following, we define the reduced control param- 
eter as 

e = r~r'f\ (5) 

Numerical solutions of Eq.(IT]) with random initial con- 
ditions for the field ?/; lead to the formation of oblique 
stripes if c > , 77 > and r > {i.e. e > \r^c''\). 
Normal stripes were never observed close to onset for the 
parameter values used in this study. 

Domain coarsening in anisotropic oblique stripe pat- 
terns was recently studied with Eq. ([T]) [s^]. This 
model has a potential structure, as it can be written as 
dtip — —5F/Sijj, with F a Lyapunov functional given by 

-2rj{d^dyyjf + c{dli:f} . (6) 

We propose here two nonrelaxational extensions of ^ 
by replacing the nonlinearity by terms that do not 
derive from a functional. We consider the cases 

NL[ij] = -V'^ - §V'(VV')^ + ^V^^VV, Model 1(7) 

Nm = £|v27A(vv')' 

Kq 

+ ^^2^(9^V)(5jV')(9«9jV'), Modem, (8) 

where ci, 02,03 are constants and = {x,y). These 

terms were previously considered in the isotropic Swift- 
Hohenberg model of Rayleigh-Benard convection [4lj (see 
also [l3|)- Expressions ([7]) and ([5]) are nonpotential when 
ci ^ — C2 and C3 ^ 2. We set ci = C2 = 1 and C3 = 3 in 
the following. 
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As shown in the next Section, all the above models 
with random initial conditions lead to the formation of 
poly crystalline configurations as the one shown in Fig. [1] 
Their large time evolution is controlled by the motion of 
grain boundaries separating domains of perfectly oriented 
zig and zag stripes. To describe patterns containing grain 
boundaries, we look for general solutions of the model 
equations of the form 29, 42] 

V'(r,i) ^A+{f,t)e'''^-^ + A-{r,t)e'''^-^ + c.c., (9) 

where — pcX ± qdj, with pc and Qc positive and given 
by Eq.([4]); c.c. means the complex conjugate. Solutions 
of the form ^ describe polycrystalline states of zig and 
zag domains with wavevectors close to the marginal ones. 
A'^ ^ and A^ ~ within a zig domain, and vice-versa. 
The amplitudes A~^ and A" are both non-vanishing in 
the vicinity of a grain boundary. 

In the limit e ^ 1, the multiple scale formalism [2^, 
l36j can be applied to the three models ([1]), ([7]) and ([5]). 
We find that, at first order in e^^^, the dynamics of yl+ 
and A^ is described by two coupled Ginzburg-Landau 
equations 

dtA+ = eA+ + ^[pldl + 2{l-Tj)paqcd,dy (10) 
+ {1 + c)qld^]A+ - 3-r\A+\^A+ - 6-/\A-\^A+ 

dtA- = eA- + ^[pldl-2{l-i^)p,q,d,dy (11) 

Kg 

+ (1 + c)qldl]A- - 37]A-]2A- - 67]^+]'^-, 

where 7 is a constant parameter whose expression de- 
pends on the model considered, 

{1, potential PK model, 

1 + ^ [tT^) ^ Model I, 
(^tStV)'' Modem. 

Equations similar to (fTU|) - PT|) were proposed as reduced 
models of oblique patterns in electroconvection [l^, |44| 
and thermoconvection [45j of nematics with planar align- 
ment close to onset. They can be recast as dtA'^ — 
-5Fgl/5A+ and dfA' = -SFgl/SA^ , where Fgl is 
a Lyapunov functional given by 

Fgl = j dr{~e{\A+\^ + \A-\^) 

4 |2 

+ ta I [Pc^x + (1 - 'n)(lcdy\A+ 1 

Kg 

4 I _|2 

+ Ta - (1 - -rDq^A \ 

Kg 

+ ^^{\dyA+? + \dyA-\^) 

+ \^{\A+\'' + \A-\^) + &-i\A+\M-?].m 



where A is the complex conjugate of A. Glose to onset 
(e ^ 1) the amplitude equations of anisotropic oblique 
striped domains are identical for the three models (fl]l. 
(O, ([H]) up to the multiplicative constant 7 (that can be 
absorbed in amplitude, time and space scales). The am- 
plitude equations have a potential structure even when 
deriving from a nonpotential model. We therefore expect 
that very close to onset oblique stripe patterns exhibit a 
generic ordering dynamics, independently of the poten- 
tial or nonpotential nature of the system. Amplitude 
equations pU|) and (fTT|) will be used to discuss some as- 
pects of the coarsening process close to onset in Section 
IV Al For a more detailed description of defect dynamics 
far from onset, higher order contributions (0(e^)) would 
be needed in the derivation of the amplitude and phase 
equations [EH^. 

III. NUMERICAL METHOD 

Equation ^ with the nonlinear terms ([7]) or ([5]) is nu- 
merically solved by using a semi-implicit pseudospectral 
method and a time integration procedure described in 
[4^ . The two-dimensional space is discretized on a square 
grid of mesh size Aa: — 1 with 1024x1024 nodes. The 
characteristic wavelength is set to Ao — 27r/fco = 8 Ax. 
The time integration scheme is stable for large values of 
the time step Ai. We set Ai = 0.5 for all the values 
of the quench depth studied. The initial condition for 
-0 is a random field with Gaussian distribution, of mean 
("0) = and variance {ijP') = e. Numerical results are av- 
eraged over 10 different initial conditions for each model 
and value of e. 

Given relation ([5]) and the fact that a random initial 
Gaussian field evolve towards oblique patterns only if 
r > 0, we choose the parameters in Eq.([T]) such that the 
smallest possible quench (e„ii„ = —ri°^) is small com- 
pared with unity. We set c = 12 and rj ~ 0.5 in all 
the calculations, corresponding to —r'c^ — 0.0196 and 
e = 11.31° (fc^/fco = 0.9901x±0.1980y). 

A typical configuration of oblique stripe domains is 
displayed in Figure [TJa) at t = 500 in the case of Model 
I for a moderate quench. As observed in electroconvec- 
tion experiments [sU . [3^ and previous numerical simula- 
tions of ^ [13] , interfaces separating zig and zag stripes 
are mainly of two kinds: horizontal chevron boundaries 
where stripes change their orientation from Q Xo —Q with- 
out any phase singularity, and slightly curved boundaries 
composed of arrays of dislocations with an orientation 
close to the vertical direction. Close to onset, isolated 
dislocations are scarcely observed within domains and do 
not seem to control the ordering dynamics. Note that iso- 
lated dislocations in anisotropic patterns are analogous 
to superfluid vortices, that are known to have a diverg- 
ing energy [13]. This is due to the fact that, unlike in 
isotropic stripes, spatial derivatives in the first order am- 
plitude equations (|10mip can be recast into a Laplacian 
form after an axis rotation. 
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(C) 



(d) 



FIG. 1: a) Oblique stripe pattern at t = 500 obtained by numerically solving Eq. (O (Model I) with e = 0.4118, c = 12, rj = 0.5 
and ci = C2 = 1; b) domain configuration of the stripe pattern in (a); c) chevron boundaries of (a); d) dislocations of (a). We 
show a fraction of 400 x 400 nodes from a total discretized system of size 1024 x 1024. 



Zig and zag domains can be identified by means of a 
(scalar) orientational local order parameter, ipd, defined 
as 



+1, if (a,v)(a,v) > (zig), 

-1, if {d,i:){dyi;) < (zag). 



The orientation field V'd of the pattern of Fig. [T^a) is 
shown in Fig. [ijb). 

Due to the anisotropy of the problem, we measure the 
growth of order along the x and y directions separately. 
For this purpose we introduce an orientational structure 
factor, defined as 



(13) 



where ^pd is the 2D Fourier transform of V-'d and the brack- 
ets indicate averages over different initial conditions. We 



define the one-dimensional structure factors along x and 
y as 



and 



Sx{qx,t) = 



Sy{qy,t) 



S{q,t)dqy, 



S{q,t)dqx, 



(14) 



(15) 



respectively. It can easily be shown that (fH|) and are 
the ID Fourier transforms of the equal time correlation 
functions along the directions x and y, respectively. 



Cxix,t)^ {ipdir ,t)^d{r +xx,t))dr, 



Cy{y,t) = / {ipd{r ,t)ipd{r +yy,t))dr 



(16) 



(17) 
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We define the correlation lengths in the x and y direc- 
tions from the inverse width of the curves {Qx ; t) a-nd 
Sy{qy,t) in Fourier space: 



Lx{t) = 



Ly{t) = 



' I^oc \'lx\Sx{qx,t)dqx' 



' iZo \ly\Sy{qy,t)dqy 
. /-!°oo'5'y(gy,t)dg^ 



(18) 



(19) 



Like in isotropic systems, the defect densities should be 
related the correlation lengths defined above if dynamical 
scaling holds. In the case of oblique stripe domains, we 
investigate separately the evolution of of chevron {pch) 
and dislocation {pdis) densities. Chevron boundaries are 
identified as the lattice nodes {x',y') where 



(20) 



The chevron density is defined as the fraction of such 
nodes, marked in black in Fig. [Ijc). We do not con- 
sider as belonging to a chevron boundary isolated hori- 
zontal clusters of nodes satisfying (|20p shorter than half 
the stripe period. Dislocations are identified by using the 
filtering procedure described in [15]. A typical configura- 
tion of dislocations is shown in Fig. [ijd), corresponding 
to the pattern of Fig. [U^a). We define numerically pdis 
as the fraction of area occupied by the black regions in 
Fig. \V[d). 



IV. RESULTS 



[ij, 122, l30|, |40| . The wavenumber k of the background 
stripes at time t is numerically obtained by averaging 
the modulus of the four wavevectors where the order pa- 
rameter structure factor {ip{q, t)ip{—q, t)) is maximum (ip 
being the Fourier transform of ip). 

Like for isotropic potential equations the selected 

wavenumber k in the PK model ([1]) is the one that min- 
imizes the Lyapunov functional After a short tran- 
sient k reaches an asymptotic value that is independent 
of e, as shown in Fig. ^ and given by the marginal 
wavenumber maximizing the growth rate ^ 



fee = \/Pc.+ <lc = \j 



2r] 



c + 2?7 — 77' 



(21) 



Nonpotential models of isotropic patterns generally ex- 
hibit nontrivial wavenumber selection far from onset, 
with k < kc- The value of k depends on e and is not 
determined by any minimization principle [3, [13, IboI . 
[isl . I49I I . We observe similar behavior for the nonpoten- 
tial anisotropic models. The selected wavenumber in the 
nonpotential models I and II is very close to kc suffi- 
ciently close to onset. As e is increased, the differences 
between the three models become more pronounced. For 
the nonpotential models I and II, k decreases with in- 
creasing e (see Fig. [2]), in qualitative agreement with 
electroconvection experiments far from onset and below 
the Lifshitz point [40|. This aspect further justifies the 
use of nonpotential equations for modeling electrocon- 
vection. The dependence of fc on e is weaker for Model I 
than for Model II, though. In Model I, k remains close to 
the marginal value (|fc — k^l/kc < 0.03) even for quench 
depths as large as e w 0.6. 



A. Wavenumber selection 



B. Correlation lengths 
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FIG. 2: Selected wavenumber as a function of e for the three 
models studied. 

Before investigating domain growth we address the is- 
sue of wavenumber selection at short times, as it can 
play an important role in the subsequent defect dynamics 



The time evolutions of the correlation lengths Lx{t) 
and Ly{t) are shown in Figs. [3] and 2] for different values 
of e and for each nonpotential model. 

Close to onset (e = 0.0392), the coarsening rates are 
well fitted by the power laws 



L,(^)^^l/^^ Ly{t)^e^'y 



(22) 



For both Models I and II, we find l/z^ = 1/zy ~ 0.45 in 
the intermediate time regime. Finite size effects become 
important at large times. These results are similar to 
those obtained for the PK potential model [l^l where 
exponents close to 1/2, the value typical of curvature 
driven growth, were reported near onset. In Section [V Al 
below, we present analytical arguments showing that the 
three models proposed can be approximatively reduced, 
close to onset only, to Model A for a non-conserved order 
parameter Q. 

At intermediate values of e, comprised at least in the 
interval [0.2, 0.4] in the case of Model I, the growth of 
both correlation lengths is well fitted by power laws with 
two distinct exponents: l/z^ ~ 0.33 and 1/zy w 0.45, as 
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FIG. 3: Correlation lengths along x and y vs. time for Model 
I at different values of e. The curves have been displaced from 
their original position for clarity. 



FIG. 4: Correlation lengths along x and y vs. time for Model 
II at different values of e. 



shown in Fig. [31 Similar results were obtained for the po- 
tential PK model [s^]. These features are also observed 
in Model II, but in a smaller interval of values of e. For 
e = 0.1176, the scahng laws ~ t°-33 and Ly - i"-^^ 
hold until t ~ 10^, see Fig. |31 This two-exponents regime 
observed at intermediate quenches in all models indicates 
that Ly grows faster than L^- However, as shown in Fig. 
[5]where both lengths are plotted together for Model I, do- 
mains are in average longer along x than along y at short 
times. Therefore, there is a characteristic time t* such 
that Lxit*) = Ly{t*). When t > t*, domains become 
longer along y than x. A qualitatively similar behavior 
is reported in electroconvection experiments performed 
in |31| , where domains are more elongated along x than 
along y at short times, although the growth along y is 
faster. 

For deep quenches, the growth rates of Lx{t) and Ly{t) 
become extremely slow in Model I, as shown in Fig. [3] 
for e = 0.8039. The transition to such arrested dynamics 
is relatively abrupt, since for e = 0.6078, Ly and Ly are 
still well described by the laws t^^^ and i^/^, respectively. 
Similar frozen states observed in the potential PK model 
were explained by the presence pinning effects acting on 



defects at large quenches [34|. However, as discussed in 
Section IV Bl freezing probably has a different origin in 
the nonpotential case. 

Whereas the ordering dynamics in Model I is similar to 
that of the potential case for a wide range of values of e, 
the nonpotential terms of Model II have dramatic effects 
at large e. For e = 0.2157 and e = 0.3137, the correlation 
lengths grow faster than a power law at intermediate and 
large times, see Fig. H) Only at short times {t < 10'^), an 
apparent scaling regime with 1/zx < 1/3 and 1/zy < 0.45 
can be observed. A secondary bifurcation from oblique 
to normal stripes actually occurs in Model II close to e « 
0.4, preventing the study of oblique domain coarsening 
for larger quenches. 



C. Defect densities 

The results presented above are further supported by 
studying the decay rate of defect densities. Assuming 
that scaling holds, the total length of chevron boundaries 
in a system of area A can be estimated as L^A/ [LxLy) ~ 
l/Ly. Similarly, the dislocation density should scale as 
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FIG. 5: Time evolution of the correlation lengths along x and 
y for Model I with e = 0.4118. The characteristic time such 
that = Ly is t* « 23000. 



l/Lx- Therefore, one should expect 

Pd^s{t) - L,{t)-^, p,h{t) - Ly{t)-^, (23) 



Pch{t) ^ t 



-1/Zah 



(24) 



with the defect exponents obeying Zdis = Zx and ZcH — 

Zy. 

Figs. [5] and [7] show the time evolution of pdis{t) and 
Pch (t) for different values of e in both nonpotential mod- 
els. 

Close to onset (e = 0.0392), the time evolution of dislo- 
cation and chevron densities is well fitted by power laws 
with 1/zdis ~ l/^c/i ~ 0.45, as shown in Fig. [5] (Model 
I) and Fig. [7] (Model II). This result further suggests 
that defect dynamics is isotropic close to onset, in agree- 
ment with the arguments exposed below in Section [V A I 
Deviations from the laws ([M)) are observed at late times 
(10'* < t) and are due to finite size effects (see below). 

By increasing e, defect dynamics becomes anisotropic 
and different laws are observed for the two types of de- 
fects. Chevron boundaries become quickly pinned and 
remain practically immobile whereas coarsening is dom- 
inated by dislocation gliding along the horizontal direc- 
tion. The reduction of the length of a chevron boundary 
can occur by the motion of the slightly curved disloca- 
tion arrays of opposite Burgers vector located on its sides 
(see Section V). This mechanism was also identified ex- 
perimentally in electroconvection [s^ ]. 

In an intermediate range of values of e, the relations 
are still fulfilled, despite that chevrons boundaries 
are pinned. As shown in Fig. [Slfor Model I, for e — 0.2157 
and e — 0.4118, Pdisit) and Pch{t) decay as power laws 
with two different exponents, z^is ~ 3 and Zch ~ 2. These 
are the same values as observed in the previous section 
for Zx and Zy, respectively. Model II leads to the same 
two exponents for e = 0.1176 and t < 10*, see Fig. [71 
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FIG. 6: Dislocation and chevron densities vs. time for Model 
I at different values of e. The curves have been displaced from 
their original position for clarity. 



Fig. [8] shows the evolution of the defect densities ob- 
tained from two system sizes, 1024 x 1024 and 512 x 512, 
with e = 0.2157. As expected, corrections to scaling 
occur sooner for the smaller system (at t w 7000), when 
the correlation lengths become comparable to the system 
size. 

Note that the decay rate for the dislocation den- 

sity agrees well with the experimental results reported in 
[32| . However, in this intermediate quench depth regime, 
our numerical results predict a much faster decay of the 
chevron density (as i^*/^) than observed experimentally 
(t^^/^ [32]). It is worth noticing that no coarsening laws 
with exponent 1/2 were ever observed in electroconvec- 
tion experiments. 

At large quenches, dislocation gliding becomes very 
slow in Model I, as shown in Fig. [S]for e — 0.8039. For 
t > 2000 the defect density saturate to an almost con- 
stant value, indicating freezing. As shown below in Scc- 
tion lVBl single defect dynamics in Model I presents some 
interesting nonrelaxational features at large quenches, 
that prevent the annihilation of small domains. 

In Model II, defect annihilation is on the contrary ac- 
celerated at large quenches and densities decay faster 



8 



10" 



10' 



"a 

Q- 10'' 



10"' 



10' 



■"□-□nm^^ fc = U.J 



£ = 0.0392 + 

£ = 0.1176 X 

£ = 0.2157 * 

£ = 0.3137 □ 



■Xx-; 




-0.28 

-0.33 
-0.33 
-0.50 



100 



1000 



10000 



100000 



10° h 



10" 



Q- 10"' 



10"- 



10"' 



□ □ 



□□[ 




£ = 0.0392 + 

£ = 0.1176 X 

£ = 0.2157 X 

£ = 0.3137 □ 



-0.50 
-0.50 
-0.50 



100 



1000 



10000 



100000 



FIG. 7: Dislocation and chevron densities vs. time for Model 
II at different values of e. 



than inverse power-laws, as shown in Fig. [7] for e = 
0.2157 and 0.3137. These results confirm those of Sec- 
tion |IVB] regarding the correlation lengths. Like for the 
other models, pinning effects might slow down defect mo- 
tion at short times (up to i ~ 1000), as dislocation and 
chevron densities decay slower than t~^^^ and t~^^^, re- 
spectively. However, nonpotential effects dominate dis- 
location motion when t > 1000 and the oblique stripe 
pattern quickly reaches fully ordered configurations. In 
Section lVBl we show that single defect dynamics is quite 
different in Models I and II far from onset. 

Note that coarsening rates faster than a i^^^ law for 
oblique stripes have been reported in electroconvection 
experiments js^]. These fast coarsening laws were ob- 
served when the background wavenumber is small, a 
property also shared by Model II (see Fig. [J). 



D. Dynamical scaling 

In the case of the potential model close to onset [13] , 
dynamical scaling properties of the structure factor were 
investigated along the directions parallel to the wavenum- 
bers and k~ . As shown above, chevron boundaries 
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FIG. 8: Late time evolution of the defect densities for Model 
I (e = 0.2157) for two different lattice sizes (1024 x 1024 and 
512 X 512 nodes). The arrows indicate the beginning of finite 
size effects. 



and dislocations have very different mobilities further 
from onset. Therefore, it sounds more natural to test one- 
dimensional scaling relations along x and y separately. 
The dynamical scaling hypothesis implies that the struc- 
ture factors (HH) and should obey the ansatz 



Sy{qy,t) = Ly{t)fy{qyLy{t)), 



(25) 
(26) 



where Lx{t) and Ly{t) arc the characteristic length scales 
previously defined through relations P^ - fTO]) : fx and fy 
are functions independent of time. 

Since the ordering dynamics of the oblique striped do- 
mains close to onset can be approximately reduced to the 
dynamics of Model A (see Section IV Ap , one should ex- 
pect that dynamical scaling holds along any direction in 
Fourier space. Fig. Oshows that the ansatz (^5)1 and ([^ 
holds over two time decades for Model II close to onset 
(e — 0.0392). The same results are observed for the PK 
potential model and Model I (not shown). However, in 
the region near fc = 0, the scaling functions in the three 
models decay rapidly (see Fig. [5]) and differ noticeably 
from that of Model A, which is quadratic for small k 
(Sll [5^. This unusual behavior may indicate that, given 
a random initial condition, oblique patterns form at early 
times different large-scale domain structures than Model 
A. 

As e is increased, isotropic dynamical scaling breaks 
down because coarsening is characterized by two grow- 
ing lengths scales with distinct exponents in the x and 
y directions. Nevertheless, Fig. [10] (top) shows that the 
one-dimensional scaling relation (|25p in the x direction 
(where Zx ~ Zdis ~ 3) holds over two time decades for 
Model II. In this regime, where similar scal- 

ing behavior is also observed along y (not shown), for 
all three models. Deviations from dynamical scaling are 
observed in the large quench depth regime, see Fig. 1101 
(bottom). 
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FIG. 10: Log-log plot of Sx{qx,t) / Lx{t) as a function of 
qxLxit) for Model II: e = 0.1176 (top) and e = 0.2157 (bot- 
tom). Solid lines correspond to the Porod's law. 



We verify that the tails of the structure factors obey 
the well known Porod's Law as illustrated in Fig. 
1101 If the one-dimensional correlation functions obey 
a scaling form at short distances, one expects that 
Sx.,viqx,y,t)/Lx,y{t) - l/{qx,yLx,y{t))''^'^^ at large qx,y, 
with d = 1 here. 



DISCUSSION 



Mapping to Model A in an external field for 
low-to-moderate quenches 



The numerical results obtained very close to onset 
show that the coarsening of oblique stripes is character- 
ized by a single exponent, Zdis ~ Zch ~ Zx ^ Zy ^ 2. This 
result can be explained by the following weakly nonlin- 
ear arguments. The amplitude equations (fTO)) and pT|) 
can be used to describe the dynamics of grain boundaries 
separating zig and zag domains close to onset (e <^ 1) in 
all models. In the absence of grain boundaries, the sta- 



tionary uniform amplitudes in a zig domain are 



(27) 



and vice- versa in a zag domain. In the vicinity of a grain 
boundary, the amplitudes change continuously from one 
uniform solution to the other along the normal coordi- 
nate. Fig. Illf a) displays the stationary numerical pro- 
files a~^{x) and a^{x) solutions of Eqs. pi]l) - pT|) close to 
onset for a vertical boundary (dislocation array). The 
sum a"*" -|- a~ and difference a"*" — a~ are also plotted in 
Fig. fllTb). We observe that a+ -I- a~ remains close to 
the constant bulk value {e/SjY^^ in the boundary region, 
while a^ — a~ changes continuously from -f (e/37)^/^ (zig- 
domain) to — (e/37)^/^ (zag-domain). For any polycrys- 
talline configuration, we define the local order parameters 



b+ = A+- 



A and cj) = — A = J ^i^d, which re- 



defines the orientational order parameter field ipd- Based 
on Fig. IllT b). we approximate </)"*" to a constant, even in 
configurations containing many interfaces: 



(28) 
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FIG. 11: a) Stationary numerical solutions {a {x),a (x)) of 
Eqs. (|10|) - (|11|) for a vertical planar grain boundary, with e = 
0.0392, c = 12, = 0.5 and 7 = 1; b) 0+ = a+ + o" and 
d>^ = a'^ — as a function of x. 



defect pinning often arises from the coupling between the 
amphtudes varying slowly in space and the short period 
modulations of the local order parameter [l^, [ssj . 
This coupling generates "non-adiabatic" terms in the am- 
plitude equations, that oscillate with the spatial coordi- 
nates and create periodic energy barriers for defects. The 
simplest generalization of Model A that incorporate sim- 
ilar pinning effects for horizontal interfaces is 



(31) 



where the periodicity of the pinning potential, of magni- 
tude p, is smaller than the interface width (27r/fcp < 1). 
The above equation describes the ordering kinetics of 
a non-conserved scalar order parameter in a stationary, 
spatially modulated external field. 

As discussed in |34l |. chevron boundaries pin at inter- 
mediate values of e while dislocations are still mobile. 
This effect can be explained by the fact that the mag- 
nitude of the pinning potential of chevrons (p) is orders 
of magnitude higher than that of dislocation arrays. A 
detailed calculation showing this and further support- 
ing the reduced model (|3ip will be presented elsewhere. 
Model A is recovered close to onset because pinning po- 
tentials tend to zero very rapidly as e — > 0, indepen- 
dently of the interface orientation [l^ \^ . Eq. ([?T|) can 
also be modified by replacing p by a term proportional 
to {ipdf'dyipd, so that the external field acts mostly on 
horizontal boundaries and not in the bulk. 



B. Nonrelaxational defect dynamics far from onset 



To find an equation for the local orientational order pa- 
rameter ^jjd, we substitute the new fields into piH) - PT|) . 
By making the change of variables 

\ c + V J V (l + c)?7 J 

(29) 

we obtain 

dr^'d = l^d + ^^^^d - ^V-d, (30) 

where T — et and — d\ + dy- This isotropic 
equation is easily recast into the well-known Model A 
for a non-conserved order parameter 0, H, 0|: pti'd — 
Tpd + V^Vd — ^d- Hence, as previously suggested [3J], the 
coarsening of oblique stripes close to onset is isotropic (af- 
ter the coordinate change (j29p ) and driven by curvature, 
with a length scale increasing as t^^^. Such regime has 
never been observed experimentally in electroconvection, 
despite that amplitude equations pU|) and pip should 
a priori describe correctly this system close enough to 
onset [H, 0,141. 

Further from onset, at intermediate values of e, the 
horizontal (chevron) grain boundaries become pinned in 
the three models considered. In pattern forming systems. 



Some qualitative insights into domain coarsening for 
Models I and II far from onset can be gained by numeri- 
cally studying the shrinking dynamics of an initially rect- 
angular zag domain embedded in a zig surrounding. The 
background wavenumber is initially fixed to the value kc- 

Fig- EUa) shows the results obtained by solving Model 
I with an embedded domain of vertical width Dy = 20Ac 
and e = 0.2157. At this moderate quench, the dynamics 
is approximately driven by the decrease of the energy 
through the reduction of total interface length. As 
mentioned in the previous Sections, chevron boundaries 
are pinned and the reduction of their length is achieved 
by the motion of the curved lateral dislocation arrays. 
When colliding, dislocations of opposite Burgers vector 
annihilate, as also observed experimentally [321] . 

Fig- HH^b) shows the evolution of the same initial do- 
main further from onset (e = 0.6078). In this case dis- 
location motion is not restricted to reduce the chevron 
length. The boundaries exhibit a more complex dynam- 
ics instead and the domain shape breaks its initial axial 
symmetry. 

To study the transition between these two dynami- 
cal regimes for Model I, we report in the upper panel 
of Fig. [T3] the evolution of the area Ad{t) of an em- 
bedded domain. At moderate values of e, Ad{t) con- 
tinuously decreases in time at a rather constant rate as 



11 




(b) 



(c) 



FIG. 12: Defect field of an initially rectangular zag domain surrounded by zig stripes, for Model I with; a) e = 0.2157 at 
t = 200,600,1800, and 3600; b) Model I with e = 0.6078 at t = 200,600,1800, and 12000; c) Model II with e = 0.2157 at 
t = 200, 600, 1800, and 2400. 



the zag domain shrinks. However, above a critical value, 
Cc ~ 0.6 for a domain of width Dy = 32Ac, Ad{t) tends 
towards a finite value at large times, suggesting that dis- 
locations effectively repel each other. The transition be- 
tween both regimes is fairly abrupt. This situation is 
reminiscent of dislocation motion modes observed exper- 
imentally, where walls of oppositely charged defects can 
approach within some distance and then move apart [3^ . 
The results above provide a picture of the arrested growth 
of order observed far from onset in Model I. The same 
single domain calculation using the PK potential model 
with e ~ 0.6078, a value close to the above Cc, shows 
that the domain completely shrinks (the bullet symbols 
of FigHni). 

We therefore conclude that, unlike in potential mod- 
els, the arrested ordering kinetics in Model I is primarily 
due to nonpotcntial repelling effects among dislocations, 
and not to dislocation pinning induced by non-adiabatic 
effects (although these might also be present). 

The lower panel of Fig. [T3] shows that, for Model I, 
the critical value Ec defined above depends on the domain 
vertical width. Interestingly, ec seems to tend towards a 



constant value (^ 0.5) as Dy —>■ oo, meaning that be- 
low this parameter value all domains shrink (coarsening 
regime). For e > 0.5, large enough domains do not shrink 
(disordered regime). This diagram suggests that Model I 
undergoes a transition to arrested coarsening as e crosses 
a critical value. 



Fig- I12f c) shows the evolution of a shrinking domain 
in Model II, for e = 0.2157. An effective dislocation re- 
pulsion is also observed, leading, in contrast to Model 
I, to the emission of pairs of oppositely charged dislo- 
cations. This mechanism gradually reduces the average 
domain width. For this value of e, the ordering kinetics 
of polycrystalline configurations is faster than a power 
law. Dislocation emission by shrinking domains at large 
quenches could be responsible for the very fast coarsening 
laws observed in Model II. 
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FIG. 13: (Color online) Upper panel: Time evolution of the 
area of an initially rectangular zag domain of fixed width 
Dy = 32Ac surrounded by zig stripes, for Model I and differ- 
ent values of e. Lower panel: Critical parameter value above 
which a domain does not shrink, as a function of the domain 
width. Below the dotted line, all domains were observed to 
shrink. 



VI. SUMMARY AND CONCLUDING 
REMARKS 

We have studied the coarsening of two-dimensional 
oblique stripe phases by using two nonpotential Swift- 
Hohenberg-like equations. The results have been com- 
pared to previous ones obtained in the potential case [13] . 

Close to onset (e ^ 1), the patterns can be described in 
all cases by two coupled amplitude equations that can be 
reduced approximately to Model A for a non-conserved 
scalar order parameter. Consequently, the coarsening 
process close to onset is self-similar and domain growth 
is characterized by a single length scale growing as t^^^ 
with z w 2. 

At intermediate values of e, a second regime is ob- 
served in both potential and nonpotential models where 



chevron boundaries become pinned and domain growth 
is driven by the horizontal motion of curved dislocation 
arrays. The one-dimensional correlation functions along 
the X and y directions still obey dynamical scaling re- 
lations, but the corresponding correlation lengths grow 
at different rates: Lj;{t) ~ t^/^" and Ly{t) ^ t^^^y, with 
Zx ^ S and Zy « 2. In this regime, the system should 
be well described by a Model A equation in an external 
modulated field (Eq.([3T|)). 

The nonpotential equations studied here reproduce 
qualitatively the observed selected wavenumbers in elec- 
troconvection experiments [i^l • The unusual but appar- 
ently robust t^^^ law for also agrees with experimental 
data [3^ . However, no 1/2 exponents were ever observed 
experimentally. 

When e ~ 0(1), a third regime can be identified, 
where defect dynamics differs noticeably in potential and 
nonpotential models. In the potential and nonpoten- 
tial I cases, coarsening rates become very slow at late 
times, the system remaining in macroscopically disor- 
dered states. Whereas freezing is due to dislocation pin- 
ning in the potential case [sj] , in Model I nonrelaxational 
effects dominate, leading to a complex dislocation dy- 
namics (quite similar to the experimental observations 
of ref. [32]) that prevents small domains to shrink. The 
study of the relaxation of single domains suggest that this 
transition to arrested coarsening in Model I is abrupt. 
This result is surprising, as it is generally a ccep ted that 
nonpotential effects enhance coarsening [illlllll. For 
Model II, nonpotential effects also affect dislocation dy- 
namics but lead to the opposite macroscopic behavior: 
defect densities decay in time faster than a power-law 
and the system orders very quickly. 

A better understanding of the transitions between the 
different coarsening regimes far from onset requires a sys- 
tematic study of the dynamics of shrinking domains and 
its role in wavenumber selection. In this spirit, a study 
of grain boundary motion in nonpotential stripe patterns 
was recently performed for isotropic systems 30]. We 
hope that the study of defect motion in nonpotential sys- 
tems far from onset will motivate new experiments. 

Finally, recent experiments ^33j studying the impact of 
noise on the coarsening rates of oblique rolls could also 
motivate future research. 



Acknowledgments 

This work was supported by CONACYT (Mexico) 
^rant 40867-F. 



[1] J. D. Gunton, M. San Miguel, and P. S. Sahni, Phase York, 1989). 

Transitions and Critical Phenomena (Academic, New [2] A. Bray, Adv. Phys. 43, 357 (1994). 



13 



P. Hohcnberg and P. Halperin, Rev. Mod. Phys. 49, 435 [29 

(1977). 

A. J. Bray and A. D. Rutenberg, Phys. Rev. E 49, R27 [30 

(1994) . 

A. J. Bray, Phys. Rev. E 58, 1508 (1998). [31 
I. M. Lifshitz and V. V. Slyozov, J. Phys. Chem. SoUds 

19, 35 (1961). [32 
S. M. Allen and J. W. Cahn, Acta Motall 27, 1085 (1979). 
M. C. Cross and P. C. Hohenberg, Rev. Mod. Phys. 65, [33 
851 (1993). 

C. Bowman and A. C. Newell, Rev. Mod. Phys. 70, 289 [34 

(1998) . [35 
J. P. GoUub and J. S. Langer, Rev. Mod. Phys. 71, S396 

(1999) . [36 
M. I. Rabinovicli, A. B. Ezorsky, and P. D. Wcidman, 

The Dynamics of Patterns (World Scientific, 2000). [37 
K. R. Elder, J. Vinals, and M. Grant, Phys. Rev. Lett. 
68, 3024 (1992). [38 
K. R. Elder, J. Vinals, and M. Grant, Phys. Rev. A 46, 
7618 (1992). [39 
M. C. Cross and D. I. Mciron, Phys. Rev. Lett. 75, 2152 

(1995) . [40 
Q. Hou, S. Sasa, and N. Goldenfeld, Physica A 239, 219 

(1997) . [41 
J. J. Christensen and A. J. Bray, Phys. Rev. E 58, 5364 

(1998) . [42' 

D. Boyer and J. Vifials, Phys. Rev. E 64, 050101(R) 

(2001) . [43 
T. Taneike, T. Nihei, and Y. Shiwa, Phys. Lett. A 303, 

212 (2002). [44 
D. Boyer and J. Vifials, Phys. Rev. E 65, 046119 (2002). [45^ 
H. Qian and G. F. Mazenko, Phys. Rev. E 67, 036102 [46 
(2003). [47 
H. Qian and G. F. Mazenko, Phys. Rev. E 68, 021109 
(2003). [48 
M. R. Paul, K-H. Chiam, M. C. Cross, and P. F. Fischer, 
Phys. Rev. Lett. 93, 064503 (2004). [49 
T. Galla and E. Moro, Phys. Rev. E 67, 035101(R) 
(2003). [50 
A. Xu, G. GonncUa, A. Lamura, G. Amati, and F. Mas- 
saioU, Europhys. Lett. 71, 651 (2005). [51 
Y. Yokojima and Y. Shiwa, Phys. Rev. E 65, 056308 

(2002) . [52 
C. Harrison et al, Science 290, 1558 (2000). [53 
C. Harrison et al, Phys. Rev. E 66, 011706 (2002). [54 
R. Ruiz, R.L. Sandstrom, and G.T. Black, Adv. Mater. 
19, 587 (2007). 



P. Manneville, Dissipative Structures and Weak Turbu- 
lence (Academic, New York, 1990). 
Z.F. Huang and J. Vinals, Phys. Rev. E 75, 056202 
(2007). 

L. Purvis and M. Dennin, Phys. Rev. Lett. 86, 5898 
(2001). 

C. Kamaga, F. Ibrahim, and M. Dennin, Phys. Rev. E 
69, 066213 (2004). 

M. Grifllth and M. Dennin, Phys. Rev. E 74, 027201 
(2006). 

D. Boyer, Phys. Rev. E 69, 066111 (2004). 

H. Q. Qian and G. F. Mazenko, Phys. Rev. E 73, 036117 
(2006). 

W. Pesch and L. Kramer, Z. Phys. B: Condens. Matter 
63, 121 (1986). 

L. Kramer and W. Pesch, Annu. Rev. Fluid. Mech., 27, 
515 (1995). 

A. Buka and L. Kramer, Pattern Formation in Liquid 
Crystals (Springer- Verlag, New York, 1996). 

M. Dennin, D. S. Gannell, and G. Ahlers, Phys. Rev. E 
57, 638 (1998). 

D. Funfschilling, B. Sammuli, and M. Dennin, Phys. Rev. 
E 67, 016207 (2003). 

H. S. Greenside and M. C. Cross, Phys. Rev. A 31, 2492 

(1985) . 

B. A. Malomed, A. A. Nepomnyashchy, and M. I. Tribel- 

sky, Phys. Rev. A 42, 7244 (1990). 

E. Bodcnschatz, W. Pesch, and L. Kramer, Physica D 
32, 135 (1988). 

M. Treiber and L. Kramer, Phys. Rev. E 58, 1973 (1998). 
E. Plant and R. Ribotta, Eur. Phys. J. B 5, 265 (1998). 

M. C. Cross, D. Meiron and Y. Tu, Chaos 4, 607 (1994). 
D. R. Nelson, Defects and Geometry tn Condensed Matter 
Physics (Cambridge University Press, Cambridge, 2002). 
A. Hari and A. A. Nepomnyashchy, Phys. Rev. E 61, 
4835 (2000). 

G. Tesauro and M. C. Cross, Phys. Rev. A 34, 1363 

(1986) . 

C. Kamaga, D. Funfschilling and M. Dennin, Phys. Rev. 
E 69, 016308 (2004). 

T. Ohta, D. Jasnow and K. Kawasaki, Phys. Rev. Lett. 
49, 1223 (1982). 

S. Puri and Y. Oono, Phys. Rev. A 38, 1542 (1988). 

Y. Pomcau, Physica D 23, 3 (1986). 

R. Gallego, M. S. Miguel, and R. Toral, Phys. Rev. E 58, 

3125 (1998). 



